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1. Introduction 

The 'MHV rules' approach proposed in Ref. [1], has led to the establishment of a 
new and powerful framework for computing large classes of previously unknown tree- 
level and one-loop scattering amplitudes in gauge theories, in a compact form, and 
without appealing to Feynman diagrams. 

In this paper, we apply the MHV rules to study the singular limits of QCD am- 
plitudes when n partons (gluons and massless quarks) are simultaneously collinear. 
This continues the program started in our earlier work [2] where MHV rules were used 
to derive multi-collinear limits of amplitudes involving only gluons. Understanding 
the infrared singular behaviour of tree-level QCD amplitudes is a prerequisite for 
computing infrared-finite cross sections at fixed order in perturbation theory. In 
general, when one or more final state particles are either soft or collinear, the am- 
plitudes factorise. The first factor in this product is a scattering amplitude that 
depends only on the remaining hard partons in the process (including any hard par- 
tons constructed from an ensemble of unresolved partons). The second factor is the 
splitting amplitude, it contains all of the singularities due to the unresolved parti- 
cles. One of the best known examples of this type of factorisation is the limit of tree 
amplitudes when two particles are collinear. This factorisation is universal and can 
be generalised to more particles [3-7] and any number of loops [8]. 

One of the main points of our approach [2] is that, in order to derive all required 
splitting functions we do not need to know the full amplitude. Out of the complete 
set of MHV-diagrams contributing to the full amplitude, only a subset will contribute 
in the multi-collinear limit. This subset includes only those MHV-diagrams where 
all of the internal propagators go on-shell in the multi-collinear limit. Moreover, 
the functions multiplying these singular propagators in the splitting amplitude are 
constrained by the MHV rules to take a purely holomorphic form: they are functions 
which depend only on the holomorphic spinor products, (i j), of the right-handed 
(undotted) spinors and not on the anti- holomorphic ones [i j]. This points towards 
a simple twistor space picture for the multi-collinear limits, in terms of a degree-one 
curve in twistor space. The MHV rules approach also enables us to calculate infinite 
sequences of splitting amplitudes - with fixed numbers of negative helicity partons 
and arbitrary numbers of positive helicity ones, or vice versa. 

The basic building blocks of the MHV rules approach [1] are the colour-ordered 
n-point vertices which are connected by scalar propagators. These MHV vertices are 
off-shell continuations of the maximally helicity-violating (MHV) n-gluon scattering 
amplitudes of Parke and Taylor [9, 10]. They contain precisely two negative helicity 
gluons. Written in terms of spinor inner products [11], they are composed entirely 
of the holomorphic products (ij), rather than their anti-holomorphic partners [ij], 
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where we introduce the common notation (piPj) = (ij) and [piPj] = [ij]. By 
connecting MHV vertices, amplitudes involving more negative helicity gluons can be 
built up. 

The MHV rules for gluons [1] have been extended to amplitudes with fermions [12]. 
New compact results for tree-level gauge-theory results for non-MHV amplitudes in- 
volving arbitrary numbers of gluons [13-15], and fermions [12, 16, 17] have been de- 
rived. They have been applied to processes involving external Higgs bosons [18, 19] 
and electroweak bosons [20]. MHV rules have also been shown to work at one-loop 
level for supersymmetric theories [24]. Building on the earlier work of Bern, Dixon, 
Dunbar and Kosower [21, 22], there has been a remarkable progress in computing 
cut-constructible multi-leg loop amplitudes in M = 4 [23-29] and M = 1 [30-34] su- 
persymmetric gauge theories. Encouraging progress has also been made using MHV 
rules for non-supersymmetric loop amplitudes [35, 36]. 

Remarkably, the expressions obtained for the infrared singular parts of Af = 4 
one-loop amplitudes (which are known to be proportional to tree-level results) were 
found to produce even more compact expressions for gluonic tree amplitudes [29, 37]. 
This observation led to the BCF recursion relations [38, 39] of Britto, Cachazo, 
Feng and Witten as well as extremely compact six-parton amplitudes [38, 40, 41]. 
These tree-level BCF recursion relations for massless particles have recently been 
generalised in two ways. In Refs. [42, 43] a new version of recursion relations was 
adopted to calculate all finite one-loop amplitudes in non-supersymmetric QCD. 
At the same time, Ref. [44] generalised BCF recursion relations to include massive 
particles at tree level. 

A comprehensive list of references and a more detailed discussion of recent devel- 
opments can be found in the recent review [45]. This progress has been stimulated 
by the original proposal of Witten in [46] of a weak-to-weak coupling duality between 
a perturbative J\f = 4 gauge theory and a topological string theory in twistor space. 

The factorisation properties of amplitudes in the infrared play several roles in 
developing higher order perturbative predictions for observable quantities. First, a 
detailed knowledge of the structure of unresolved emission enables phase space inte- 
grations to be organised such that the infrared singularities due to soft or collinear 
emission can be analytically subtracted at NLO [48-50] or at NNLO [51]. Second, 
they enable large logarithmic corrections to be identified and resummed. Third, the 
collinear limit plays a crucial role in the unitarity-based method for loop calcula- 
tions [21, 22, 52, 53]. 

In general, to compute a cross section at N n LO, one requires detailed knowledge 
of the infrared factorisation functions describing the unresolved configurations for 
n-particles at tree-level, (n— l)-particles at one-loop etc. The universal behaviour in 
the double collinear limit is well known at tree- level (see for example Refs. [54, 55]), 
one-loop [21, 56-60] and at two-loops [61, 62]. Similarly, the triple collinear limit has 
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been studied at tree-level [3-6] and, in the case of distinct quarks, at one-loop [63]. 
Finally, the tree-level quadruple gluon collinear limit was derived in Ref. [2, 7]. 

Our paper is organised as follows. In Section 2, we briefly review the colour 
ordered formalism that underpins the MHV rules. The relevant MHV vertices are 
given in Section 3. Section 4 describes the procedure for taking the collinear limit 
while the analytic structure of the splitting functions is discussed in Section 5. We 
write down general collinear factorization formulae in Section 6, which are valid for 
specific numbers of negative helicity partons and an arbitrary number of positive 
helicity partons. These results involve quarks and gluons in the collinear set and 
are complementary to the multi-gluon splitting functions derived in Ref. [2]. Specific 
explicit results for the collinear limits of up to three collinear partons are given in 
Sec. 7. Our findings are summarized in Sec. 8. 

2. Colour-ordered amplitudes 

Tree-level multi-particle amplitudes can be decomposed into colour-ordered partial 
amplitudes. For gluons only, this decomposition is given by 

A n {{p t ,X l ,a l })=ig n ' 2 Tr(TMD...TM.))A n (a(l Al ,...,n A ")). (2.1) 

Here S n /Z n is the group of non-cyclic permutations on n symbols, and j Xj labels 
the momentum pj and helicity A-,- of the j th gluon, which carries the adjoint repre- 
sentation index aj. The T ai are fundamental representation SU(iV c ) colour matrices, 
normalized so that Tr(T a T b ) = 5 ab . The strong coupling constant is a s = g 2 /(47r). 
Note that the MHV rules method of Ref. [1] is used to evaluate only the purely 
kinematic amplitudes A n . Full amplitudes are then determined uniquely from the 
kinematic part A n , and the known expressions for the colour traces. 

For processes involving a quark- ant iquark pair and an arbitrary number of glu- 
ons, the colour decomposition is given by 

A n ({pi, Xi, di}, {pj, Xj, ij}) (2.2) 

= ig n ~ 2 (™ 2) • • •™- 1) )uin A n (l x q \a(2 x \ . . . , {n - l)*- 1 ),^") > 

<TgS„_2 

where S , „_ 2 is the set of permutations of (n — 2) gluons and the fermions carry 
the fundamental colour labels %\ and i n . By current conservation, the quark and 
antiquark helicities are related such that Ai = — A n = A where A = ±|. 

When an additional photon with momentum P 1 is emitted, the amplitudes have 
the following form, 

A n ({pi, Xi, a^, {pj, Xj, ij}, P T ) (2.3) 
= ieg n ~ 2 (™ 2) • ■■T a ^) nin A n (l x q \a(2 x \ ...,(n- l)^" 1 ), n A "; P 7 ) , 
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where e is the electric charge of the quark. 

When there are two quark- ant iquark pairs the tree-level amplitude can be de- 
composed into colour ordered amplitudes as, 



n-4 , 

A n ({pi, \i, a*}, {pj, Xj, ij}) = ig n ~ 2 X] | 

(TMD . . ■T a ^) ilin (T a pW . --T a ^) ls+lls 

xA n (l x , (7(1), . . . , a(k)), s^ x '; (s + 1)£, p(l), p(l),nf) 

xA n (l x , (7(1), . . . , a(k), sf; (s + 1)* p(l), . . . , p(l), nf) } (2.4) 

where S k and Si are permutation groups such that k + I = n — 4 and represent the 
possible ways of distributing the gluons in a colour ordered way between the quarks. 
For i — j — 0, (T fll . . .T a i) k i reduces to 5 k i- We see that the two amplitudes A n 
and A n correspond to different ways of connecting the fundamental colour charges. 
For the A amplitudes, there is a colour line connecting q and Q and a second line 
connecting Q and q, while for the QED-like A amplitudes the colour lines connect q 
to q and Q to Q. Any number of gluons may be radiated from each colour line. As 
before, by current conservation, the quark and antiquark helicities are related such 
that X q = —Xg = A and Xq = —Xq = X' where A, A' = ±|. 

3. MHV amplitudes 

The colour ordered n-gluon MHV amplitude is given by 

A n (l + , . . . , m 1 , . . . , m 2 i ■ ■ ■ , n + ) = =Jj / • • , i \ ' (3-1) 
while the two-quark multi-gluon MHV amplitudes are, 

/ -i \2-2A / \2+2A 

a /- \ \^ \ m 1) (run) 

^ '"-■■■-"« 7 >- TE.J, + i) ■ (3 ' 2) 

Here the helicity of the quark is denoted by A = ±| while . . . denotes an arbitrary 
number of positive helicity gluons. Amplitudes for a quark-antiquark pair, many 
gluons and a photon are given by, 

A(l x ,...,n q x ;P J ) = (pi)(12) ... (nP) . ( 3 - 3 ) 

/ -i \ 2— 2A / \ 2+2A 

. . . ,m , . . . P+) = L-L_^-i_ . (3.4) 
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In the four-quark case, there are four MHV amplitudes where two of the fermions 
have negative helicity and two have positive helicity for each colour structure. For 
each helicity configuration we can write, 

, • • • , Sq, (s + 1)q, ...,riq) = - j=p |^ ^_|_ ^ "> (3-5) 



Mit ■ 1)5. ■ ■ ■ > ) = nl^ + i) ' (3 ' 6) 

A(l- sUs + D- n+) <^)<^ + 1)" + (37) 

^nl-Lg , ■ ■ ■ ,Sq,{S -\- L)q, . . . ,fl q ) — YT (ll + l) ' 



(1 s) 3 (ns + 1) 

A n (l", ...,8q,{ 8 + 1) J, ..-,<)= ^^^^^ (3-8) 
with the other colour ordering given by, 



-*4ji(1^~> • • • ) S q 1 ( S + 1)q> • • • ) n Q) — j-jn ^^_|_^ ' (3-9) 



(In) (ss+ l) 3 



..., S r,( s+ l)-,..., n +) = V\„ ; / v , (3.10) 



AnU g , • • • ,S q , {S+ L) Q , . . . ,TIq) - Yl U + ' 

i n (l", . . . , 4, ( a + 1)5, . . . , n~ Q ) = (3-12) 



The MHV amplitudes are related by parity and can be obtained by conjugating 
the MHV expressions, 

A n (l x \ n A ") = (-1)" {A n (r x \ n"*"))* , (3.13) 

and similarly for the A amplitudes. 



4. Collinear limits 

To find the splitting functions we work with the colour stripped amplitudes. For 
these colour ordered amplitudes, it is known that when the collinear particles are not 
adjacent there is no collinear divergence [7]. Therefore, without loss of generality, we 
can take particles 1 . . . n collinear. 

The multiple collinear limit is approached when the momenta pi, ■ ■ ■ ,p n become 
parallel. This implies that all the particle subenergies s^- = (pi +pj) 2 , with i,j = 
1, ... are simultaneously small. We thus introduce a pair of light-like momenta 
P u and f " (P 2 = 0, £ 2 = 0), and we write 

(Pl + ---+Pny = P U + jf^p, S i , j = (p i + ..-+ Pj ) 2 , (4.1) 
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where si jW is the total invariant mass of the system of collinear partons. In the 
collinear limit, the vector P v denotes the collinear direction, and the individual 
collinear momenta are "p\ — > ZiP v . Here the longitudinal-momentum fractions Zi are 
given by 



(4.2) 



and fulfil the constraint YlT=i Zi = 1- To be definite, in the rest of the paper we work 
in the time-like region so that (s^- > 0, 1 > Zi > 0). 




N 




Figure 1: Factorisation of an iV-point colour ordered amplitude with gluons p±, . . . ,p„ 



collinear into splitting function for P 
amplitude. 



1, . . . ,n multiplied by an (N — n + l)-point 



As illustrated in Fig. 1, in the multi-collinear limit an iV-particle colour ordered 
tree amplitude factorises and can be written as 

A N (1 X \ N x ") - split(l Al , . . . , n A " - P x ) x A N _ n+1 ((n + 1) A ^, . . . , N x » , P x ). 

(4.3) 

This labelling of the splitting amplitude split(l Al , . . . , n Xn — > P x ) differs from the 
usual definition because we use the momentum and helicity that participates in the 
resultant amplitude P x rather than —P~ x . With this choice, it is easier to see how 
the helicity is conserved in the splitting, i.e. helicity A 1 ,..., A" is replaced by A. 
Since eq. (4.3) applies for all N, we can use it to derive the splitting amplitude by 
systematically choosing iV = 3 + n. In this case, we always factorise onto a four-point 
amplitude. 



5. Analytic structure of splitting amplitudes 

The MHV rules of Ref. [1] were developed for calculating purely gluonic amplitudes 
at tree level and extended to amplitudes involving fermions in Ref. [12]. In this 
approach all non-MHV iV-particle amplitudes (including MHV) are expressed as 
sums of tree diagrams in an effective scalar perturbation theory. The vertices in this 
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theory are the MHV amplitudes of Eq. (1.1) continued off-shell and connected by 
scalar propagators 1/q 2 . 

Following [2], we classify collinear limits according to the difference between the 
number of negative helicity particles before taking the collinear limit, and the num- 
ber after, AM. Splitting amplitudes are calculated using the factorisation formula 
eq. (4.3). To facilitate the calculation, it makes sense to factorise onto hard ampli- 
tudes with the simplest analytic structure. Hence, in the MHV-rules formalism we 
will always factorise onto MHV amplitudes which are listed in section 3. In this case 
we find that AM of the splitting amplitude satisfies the relation, 



AM + 2 = N. 



(5.1) 



where 2 is the number of negative helicities in the hard MHV amplitude, and N_ 
is the total number of negative helicities in the full amplitude. AM determines the 
order of MHV diagram [1] for the full amplitude An 
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A N = MHV 



A N = NMHV 



A N = NNMHV 



(5.2) 



and so on for all AM > 2 cases. 



If we choose to use MHV rules, we extract the splitting function by factorising 
onto MHV amplitudes. Splitting amplitudes are then classified by the difference in 
the number of positive helicity particles, AP, and similar observations apply. 

In general, any splitting amplitude can be obtained from either MHV or MHV 
rules. A simple power counting argument [2] gives 

1 



split oc 



j AM / ^AP 



(5.3) 



For an MHV-rules diagram to contribute to AM ^ collinear limits, it must contain 
anti-holomorphic spinor products [ij] of collinear momenta. However, because on- 
shell MHV vertices are entirely holomorphic, within the MHV rules there are only 
two potential sources of the anti-holomorphic spinor products. One source is scalar 
propagators l/s^ = 1/ (ij) [j i] which connect MHV vertices. The second source is 
the off-shell continuation of the corresponding connected legs in the MHV vertices. 
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Each off-shell continued leg of momentum P gives rise to a factor (iP) oc (i|-P|?7] 
which amounts to anti-holomorphic factors of the form [jrf\. When the reference 
spinors rja are kept general, the 77-dependence must cancel and therefore the off-shell 
continuation cannot give rise to an overall factor of [ij]. 

This implies that within the MHV rules, the anti-holomorphic spinor products 
in (5.3) arise solely from the internal propagators. Since AM = vmhv ~ I, 1 where 
vmhv is the number of MHV vertices in the diagram, the total number of internal 
propagators is AM, in agreement with (5.3). Similarly, in the MHV approach, 
the holomorphic products would arise solely from internal propagators whose total 
number in MHV diagrams is AP. 

More precisely, it follows that all splitting amplitudes can be recast as 

s P lit = EfOT — /«» ( 5 - 4 ) 

1 U,j=l Si ,3 

= E /(ID ( 5 - 5 ) 

1 U,j=i s * ,3 

where the first expression follows from the MHV rules representation, and the second 
expression - from the MHV formalism. Here the summations are over all inequiv- 
alent choices of AM (AP) products of vanishing kinematic invariants s^j which 
corresponds to different MHV (MHV) rules diagrams. The coefficient functions / 
depend only on holomorphic spinor products, while the MHV coefficients / are purely 
anti-holomorphic. Moreover, / and / have dimensions, 

1 1 

f ^ ^ ^ AP— AM ' / ^ j j AM— AP • (^) 

The fact that / (/) is purely (anti)-holomorphic suggests a simple twistor-space inter- 
pretation. All splitting functions can be represented as sums over the corresponding 
poles in s with the coefficients being supported on a single degree-one curve in (anti)- 
twistor space. This pure (anti)-holomorphic representation of multi-collinear limits 
is specific to the MHV (MHV) formalism and is lost in the usual Feynman- diagram- 
type approaches as in Ref. [7], or in the BCF recursive approach, as shown in [2]. 

We further note that MHV rules for collinear limits are substantially simpler 
than the rules for the full amplitudes. Collinear splitting functions follow from a 
subset of the MHV rules diagrams [2]. The subset is determined by requiring that 
all internal propagators are on-shell in the multi-collinear limit. 2 This is a powerful 
constraint on the types of the contributing diagrams and it simplifies taking the 
collinear limit dramatically. 

In principle, AM = Vmhv ~ v 'mhv wnere v 'mhv ^ s ^ ne number of MHV vertices remaining in the 
factored amplitude. The splitting function is independent of v' MHV , and because we systematically 
choose to factor directly onto a single MHV vertex, we set v' MHV = 1. 

2 This is dictated by cq. (5.4) in the MHV formalism (or eq. (5.5) for MHV rules). 
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As mentioned earlier, each splitting amplitude can be calculated in both the 
MHV and in the MHV approaches. In practice, eqs. (5.3)-(5.5) imply that the MHV 
approach is simpler if AM < AP, while the MHV approach is more compact in the 
opposite case, AP < AM. 

In most of what follows we will concentrate on the splitting amplitudes with 
AM < AP and will follow the MHV rules. The remaining amplitudes with AP < 
AM are obtained from these by complex conjugation. 

5.1 An example 



When AM = AP both MHV and MHV rules are expected to yield results of similar 
complexity. As an example, let us consider a triple collinear splitting with AM = 
AP =1. In full generality, the MHV (MHV) rules approach should generate a 
maximum of three terms corresponding to simple poles in Si^, £2,3 and 51,3 = (pi + 
P2+P3) 2 - For the specific splitting 1~, 2^, 3 Q — > P~, the MHV rules approach yields, 

(2 3) z 2 



split(l+ 2+3^P+) 



while the MHV rules approach finds, 
split(l+ 2+ 3q->R h 



S2,3«12) v ^+(13) v ^)(2 2 + Z 3 ) 

((13)^+ (2 3) Jztf 
"s li3 ((12) v ^+(13) v ^) (23) ' 

ziz 3 [2 3] 



(5.7) 



[12f 

Sl)3 ([12] v ^+[13] v ^)[2 3] ' 



(5.f 



As expected, the s 12 pole is absent because there is no qQ collinear limit. By taking 
the limit of a Feynman diagram calculation, Ref. [7] finds, 

splits, 2J3 5 P?) = + ^MiW3 . (, 9) 

q Q V 82,3(^2 + ^3) Si i3 S 2 ,3 

Results (5.7), (5.8) and (5.9) are for the same amplitude and all three expressions 
agree numerically. But the analytic form of these specific representations is different. 
In agreement with eqs. (5.4)-(5.5), the functions accompanying the 1/s poles, are 
holomorphic in the MHV result (5.7), are anti-holomorphic in the MHV expression 
(5.8), while the Feynman diagram result (5.9) contains a mixture of holomorphic and 
anti-holomorphic terms. (In this case, it happens to give a more compact result.) In 
general, the limit of an amplitude computed using the BCF recursion relations will 
also provide a mixed holomorphic/anti-holomorphic splitting function (as discussed 
in sections 4.2.2 and 4.2.3 of Ref. [2]). In this specific case, taking the collinear 
limit of the compact expression for the appropriate six-parton amplitude given in 
Ref. [41] exactly reproduces the MHV result of eq. (5.8). 
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Figure 2: MHV topologies contributing to (a) Split + (mi) and (b) Split_ (mi, 771,2). Neg- 
ative helicity particles are indicated by solid lines, while arbitrary numbers of positive 
helicity particles emitted from each vertex are shown as dotted arcs. All particles that are 
not in the collinear set must be emitted from the left-hand vertex. 

6. General results 

In this section we give the results for the multiple collinear limit of quarks and gluons. 
We categorise the results according to the number of quarks involved in the limit. 
In each case, we give the general results for collinear limits with AM =0,1 and 
involving an arbitrary number of positive helicity particles. 

Limits of the type split(l + , . . . , n + — > P + ) and split (1~, 2 + , . . . , n + — > P~) can 
contribute to the AM = 0, and these collinear splitting functions are straightforward 
to derive directly from the simple MHV vertex. 

For the remaining splitting functions, it is useful to introduce the more compact 
notation 

split(l + , . . . , rrii , . . . , mj, • • • , m~ , . . . , n + — > P ± ) = Split ± (m 1 , . . . , m r ) . (6.1) 

For AM = 1, there are two possible types of splitting function, Split + (mi) and 
Split _ (mi, . . . , m r ). The possible MHV topologies contributing to these splitting 
functions are illustrated in Fig. 2. Only negative helicity particles are shown. In 
the collinear limit, the propagator goes on-shell. Any MHV diagram with a hard 
particles emitted from both vertices produces an off-shell propagator. This means 
that only particles from the collinear set are allowed to couple to the right-hand 
vertex. All hard partons couple to the left-hand vertex. 

Throughout we adopt the notation of Ref. [2]. In order that the limits can be 
read directly from the MHV diagrams, we make the following substitutions. If a is 
a particle from the collinear set, b is a particle which is not in the collinear set, and 
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q is the sum of the collinear momenta from i + 1 to j, then 

j 

(aq) - [Pv] E ( al ) v^i = l P v] A(i, j; a), (6.2) 

J=i+1 

(6g> - [Pi7](6P> £ *i. ( 6 - 3 ) 
i=i+i 

(ba) - (bP)y/^. (6.4) 
The A is defined as 

j 

A(t,j;a) = £ (aOv^, (6-5) 

J=i+1 

noting that the boundary terms involving either (0 1) or (nn + 1), are given by, 



nn + 1) ^ (6.6) 



A(i,n;n + 1) EL+i^' 
(01) ^ 



We also introduce 



(6.7) 



,2 



D(i,j,q i+ld ) = q ^' j A(i,j;i)A(i,j;i + l)A(i,j;j)A(i,j;j + l) . 

(6.8) 

6.1 One quark in the collinear set: q(ng) — > g 
6.1.1 AM = 

This is the simplest case which is read directly off the single MHV vertex. For 
positive helicity quarks, we use the two-quark MHV amplitude of Eq. (3.2) and find, 



split(l+, . . . , n+ - P+) = V* — ■ (6.9) 

For negative helicity quarks, 

spht(l-, . . . , n+ - P~) = . (6.10) 

Note that helicity conservation ensures that the helicity of P is the same as that of 
q. It is often convenient to combine results for quarks of helicity A = ±| such that, 

r— 2-2A 

split(l^,...,n + ^P A ) = ¥J_ . (6.11) 

r \ q i i q) . i-rn— 1 /; ; i -I \ \ / 
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Using parity we find, 

split (1+ ...,n 



(-i) 



n-1 



■2+2A 



(6.12) 



\/ziz n n^i 1 [i, i + 1] 

The amplitudes where an antiquark is collinear with several gluons are obtained 
by charge conjugation. 

6.1.2 AM = 1 

Because of helicity conservation, AM = 1 implies that a single gluon has negative 
helicity. When the quark has positive helicity, then the MHV diagrams contributing 
in the collinear limit correspond to topology (a) of Fig. 2. There are two types of 
diagram - one class where the quark is emitted from the right-hand vertex (and the 
propagating particle is a quark) and one class mediated by gluon exchange where 
the quark is emitted from the left-hand vertex. We find, 

,n + ->P+) ' 



s P iit(r 



m 



x 



A 3 (0,j;m)(lm> 



m— 1 n 



1=1 jr=m 



D(iJ,q i+hj ) 



(6.13) 



In the same manner, for negative helicity quarks, the allowed MHV diagrams corre- 
spond to the first and second topologies shown in Fig. 2(b), 

1 



split (1" 



m 



x 



+ 



-E 



/ „• \ 3 

A(0, j; m)(l m) 3 



m— 1 n 



i=l j=m 



A 4 (i, j;m) 
D(i,j,q i+ i,j) 



(6.14) 



6.2 Two quarks in the collinear set: (ng)qq — > g 

In this collinear limit, the gg pair is in the adjoint representation and effectively acts 
as a gluon. 

6.2.1 AM = 

This is the simplest case which is read directly off the single MHV vertex. Unlike the 
previous case, here we start with a two-quark MHV amplitude and factorise onto a 
gluonic MHV amplitude. Alternatively, we could start with a four-quark amplitude 
and factorise onto a two-quark amplitude. For quarks with helicity A = ±|, we find, 



spUt(l+...,s^( s +l)*.,n + -►/>") 



■2+2A 



-2-2A 



1 



(6.15) 
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6.2.2 AM = 1 

For amplitudes of the Split + (mi)-type, we find 



split(l+,..., S r\( s + l)J, 



,n 

s— 1 n 



j2 j2 A2+2A (^^ g ) A2 ~ 2A (^^ s + 1 ) (G 1C) 



j=s+l 



There are four diagrams contributing to splitting functions of Split_(mi, 7712) 
type 3 , 



split(l + ,...,sr A ,( s +l)J,...,m-,...,n + ) = 



1 



sjz^n nr=i 1 1 + i ) 



X 



s— 1 n 

EE 

i=0 j=m 



(sm) 2+2A (s + lm) 
D(i,j,q i+1J ) 



,2-2A / 



Zfc 

vfc=i+l 

+ g g A 2+2A (7,jj g )A 2 - 2A (7,j; S + l) c 2 

i=0 j'=s+l 



D(i,j,q i+ i,j) 



-E 

j=m 
m— 1 n 

+ E E 

j=s+l j=m 



A 2+2A (s, j;m)(s + lm) 



2-2A 



2-2A 



-2+A 



A 4 (*,j;m) 
D(i,j,q i+ i,j) 



E 

v fc=s+l 



2+2A 



\/ z s+l 



2-2A 



(6.17) 



Splitting functions of the type split(l + , . . . , m~, . . . , Sq ,(s+ 1) A , . . . , n + — > P~) are 
obtained by line reversal. These results for the two-quark sector are sufficient to 
calculate all splitting amplitudes for up to four partons. 

6.3 Two quarks in the collinear set: q(ng)q — > 7 

In this collinear limit, the q . . .q system forms a colour singlet and effectively acts as 
a photon. 



6.3.1 AM = 

In this limit the four-quark A MHV amplitudes of eqs. (3.9)-(3.12) factorise directly 
onto the two-quark+photon amplitudes of (3.3). We find that, 



split(l A ..,n, 



(12) ■■■(n- In)' 



(6.18) 



3 Diagrams where both the negative helicity fermion and gluon couple to the right-hand vertex 
in Fig. 2(b) can be mediated by either fermion or gluon exchange. 
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-f •" 

N 
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m 
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71 + 1 



AT 



n 



m S m 

n + 1 n 



Figure 3: MHV topologies contributing to the two quark collinear limit of the type 



split (1; 



, m 



P ). Quarks of type Q (q) are shown as green(red)-dotdashed 



lines and negative helicity gluons as black solid lines. The negative helicity photon is shown 
blue dashed line. 



6.3.2 AM = 1 

For amplitudes of the Split, (mi) type, there is a single MHV diagram and we find 



split (1 ,...,ni 



V^Uj=i (h 1 + 1) 



A(0,n; l) 3 A(0,n;n) 
D(0, n; gi n ) 



(6.19) 



As in the previous case, there are four diagrams shown in Fig. 3 contributing to 
splitting functions of Split _ (mi, m 2 ) type such that, 



split (1 ,. ..,m ,...,nt -> P- 



v z ^ z n nr=i 1 ( z > 1 + i ) 



X 



r TO— 1 n— 1 

EE 



TO— 1 



i=l j=m 

(lm)" 1 A(0, j; m) 



\ (rim) A(z, n; m) 3 3 

^ D(i,n;qi +ltn ) Zl 



E 



nlJ^Zk) + 



(1 m) 3 (nm) 
D(0, n; gi, n ) 



v /c=i+l 

(6.20) 



6.4 Three quarks in the collinear set: q(ng)QQ — > g 

In this configuration, the Q is adjacent with Q and therefore the vertices in the MHV 
rules include the four-quark amplitudes of eqs. (3.5)-(3.8). The factorised amplitude 
is a two-quark MHV as given in eq. (3.2) Furthermore, since the helicity of q is 
conserved and the helicities of Q and Q are opposite, there are no AM = splitting 
functions. 
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6.4.1 AM = 1 

For AM = 1, the two diagrams (with quark and gluons exchanged) of Split + (rai)- 
type yield, 



split(l+, . . . , s Q , (s + 1) J, . . . , 



n + -> P+) = 



x 



Pi ^ A(i,j;s + l)A 3 (i,j;s) 



j=s+l 

split(l+,...,sj,(s + l)g,...,n 



_ ^ (l g )A 2 (0,j; S )A(0,j; g +l) 



a=l 



X 



V Z l Z n Wl=l (I, 1 + 1) 

yi y^ A 3 (i,j;s + l)A(i, j; j 



^ (l g )A 3 (0,j; g +l) \ 

~h \h Zk ) 

Similarly, the two diagrams of Split_(mi, m 2 )-type yield, 

i 

split(l~, ...,Sq,(s+ 1) J, . . . , n + -> P") = 



\f z ^YCi=l (i, 1 + 1) 

>p A(i,j; s + l)A 3 (i, j;s) 3 
_ ^, (ls) 3 A(0,j;s+l) V 



split (l g ,...,sj, (s + l) Q , ...,n 

s— 1 rt 



P7) = 



x 



V z i z ^Y[ r i=i (h 1 + 1) 

A 3 (i,j;g + l)A(i,j;g) 3 



(6.21) 



(6.22) 



(6.23) 



_ y. (l S )(l g +l) 2 A(0,j; S +l) 



j=s+l 



Y. Zk 



.(6.24) 



6.5 Three quarks in the collinear set: q(ng)qQ — > Q 

Here the relevant vertices in the MHV rules include the A four-quark amplitudes 
of eqs. (3.9)-(3.12) and the factorised amplitude is a two-quark MHV as given in 
eq. (3.2) As in the previous case, the quark helicities are constrained such that there 
are no AM = splitting functions. 
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6.5.1 AM = 1 

There are two diagrams for both Split + (mi)- and Split _ (mi, m 2 )-types and we find, 

1 



split(l+,..., S r,( s + l)+,...,n+^P+) = 



x 



A(0,s;l)A 3 (0, S ; S ) 



D(0,s,q hs ) 



-\J z s+\ 



^ ( gS + l)A 2 (0,j; S )A(0,j;l) (1 
^ D(0,j,q ltj ) \^ Zk 



j=s+l 



,fc=l 



, (6.25) 



split(l,,...,s+ (s + l)+ ...,n + ^P+) 



x 



A 3 (0,s; l)A(0,s;s) 



£>(0, s,gi, s ) 



^ ( Sg + l)A 3 (0,j;l) / * 



j=s+l 



,fc=l 



(6.26) 



and, 



splitd:, ...,.-„(.+ 1) Q . ...,„♦-. F Q ) = ^_^ {lJ + i) 



X 



A 3 (0, S ;l)A(0, S ; g ) 3 

D(0,s, qi , s ) VZs+1 



D(0,j,q ld 

split(l~ . . . , s+ (s + 1)5, . . . , n + -> Pq) 



,fc=i 
1 



(6.27) 



J'=s+1 



X 



A 3 (0, S ;1)A(0, S;S ) 3 

^(0, S ,g M ) 



A ( ss + l)(l s +l) 2 A(0,j;l) , 
^ D(0,j, qi>j ) \^ 



j=s+l 



Zk 



.(6.28) 



6.6 Four quarks in the collinear set: QQ(ng)qq — > <? 



This limit is associated with the four-quark A-type colour ordered amplitude and is 
obtained by factoring onto a gluonic MHV. 



6.6.1 AM = 1 

Because of helicity conservation for the quarks, AM = is forbidden. Furthermore, 
at least two negative helicity quarks participate in the scattering so that AM = 1 
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Figure 4: MHV topologies contributing to the four quark collinear limit of the type 
split(l+,...,s£,(s + + l)- x ',...,n+ -> P-). Quarks of type Q (q) are 

shown as green(red)-dotdashed lines and negative helicity gluons as black solid lines. 



splittings must be of the Split_(m 1; m 2 )-type. The five contributing diagrams are 
shown in Fig. 4. Explicit evaluation of the four independent helicity configurations 
yields, 



split (1" 



split (r 



q 5 • • • 5 



v z i z n nr=i 1 * + !) 



-i t-i 



A(«,j;s)A 3 (i,j;s + 1) 3 

2 , 2 / nZ- a ~ — ^ v Zt V z w 



i=0 j=s+l 
t-l n 



D{i,j,q i+ i,j) 



i=s+lj=t+l 
s— 1 rt 



D{i,j,q i+ i,j) 



j=0 j=t+ 



x D(i,j,q i+1J ) 



K k=i+1 



^ A(g,j;t + l)( g + lt) 3 / ^ 



s-1 



i=0 



D(s,j,q s+1;j ] 

A(i,t;s){s + lt} 3 
D(i,t,q i+ i jt ) 



\k=s+l 

' :, - J 1 Zfc 
v fc=i+l 



(6.29) 



-1 t-l 



A 3 (z,j;s)A(i,j;s + 1) _„ _ 

Z^ ruT~rz — \ V^t+l 



■3 



i=0 j=s+l 



D(i,j,q i+ i,j) 



i=s+l j=t+l 



D(i,j,q i+ i,j) 



- 17- 



i=Oj=t+i "\i,j,q t+li3 ) \ k=i+1 / 
^ A(s,j;t+l)A*(s,j;t)(s + lt) 3 ( 



z k 



j=t+i v 
^ A{t,t;s)(s + lt)(ts) 2 f A 



3 



Zk 



(6.30) 



S P lit(l + ,...,4,( S +l)Q,...t + 

i=S+lj=t+l V U)>iltlj/ 



;=s+i j=t+l 

*p A (s+lt)(gt+l)( g +lt+l) 2 / ^ 



+EE 

i=0 j=t+l 



D(i,j,q i+ i,j) 

">V lt+1> ^ f E „ N 



split (1+,. 



_ A( S ,j;t+l)( S +^ / ^ 
D(s,j,q s+1J ) 

^ A(t,t;s)A 2 (t,t;s + l)(s + lt) 3 / * 

i 



i=0 



_ r> 1 + 1 



(6.31) 



g e A3( ^ )A( ^r 1} ^ 

t^-fr^ D(i,j,q i+1J ) 



^ ^ A(z,j;t)A 3 (i,j;t+l) , 

i= 0j =t+i "Kh^Qi+ij) \ k=i+1 J 
" A 3 (., It) , 3 /U, 



(6.32) 
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Figure 5: MHV topologies contributing to the four quark collinear limit of the type 
split(l+,... ,s|,(s + 1)- A ',...i£',(t + . . . , ra+ -> P~). Quarks of type Q (q) are 
shown as green(red)-dotdashed lines and negative helicity gluons as black solid lines. 



6.7 Four quarks in the collinear set: Qq(ng)qQ — ► g 

This limit is associated with the four-quark A-type colour ordered amplitude and is 
obtained by factoring onto a gluonic MHV. 

6.7.1 AM = 1 

As in the previous case, helicity conservation for the quarks, ensures that AM = 
is forbidden and that Split + (mi)-type AM = 1 splittings are absent. The four 
contributing diagrams of Split_(mi, 777.2) -type are shown in Fig. 5. 
The four independent helicity configurations are given by, 

.,/.(/ + nr, >r - . 1 



split (r 



tg,{t + l 



Q) • • • J 

s— 1 n 



sjzizn nr=i 1 {h 1 + 



E E 

L j=0 J=t+1 



(t + lt) 3 (gg + l) 

D(i,j,q i+ i,j) 



Zk 



^ A(s,j; g + l)(t + lt) 3 

- 2^ — ; — v: 



j=t+i 

s-1 



D[s,j,q s+ i 



z * 

yk = S + l 



^ A 3 (z,t;t)( gg + 1) E 

^ D(z, £, g i+ i )t ) 

A 3 (s, t;t) A(s,t;s + 1) 



D(s,t,q s+h 



E 

fc=i+l 
3 



2fc 



(6.33) 
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' ,/ ) V- ' ^Q> • • • 5 

s— 1 n 



n 



p- 



\/ZlZn YYl=l ( l > 1 + l ) 



X 



EE 

i=a j=t+i 



(t + lt)(ss + l){s + lt+l)< 
D(i,j,q i+ i,j) 



Zk 

y k=i+l 
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- y A( g ,j; g +l)(t + lt)( g + lt+l) 2 / ' V 



s-l 



^ A(i,t;t)A 2 (i,t; s + l)(ss + 1) 



i=0 



D(i,t,q i+ht ) 



fc=z+l 



split(l + , . . . , Sq, (s + 1)+ 



A(s,t;t)A 3 (s,t;s + 1) 
D(s, f, q s+ i, t ) 

..tr,(t + l)+,...,n+^p-) 



(6.34) 



v^i^ nr=i 1 * + 



s-l 



E E 

i=0 j=t+l 



(t+ lt)(ss + l)(ts) 2 ( z 



\k=i+l 



_ y A(s,j;s+l)A 2 (s,j;t)(t+lt} 
- ! A(^t;t)( Sg + l)(t g ) 2 /—/^ . 

_ . y/Zt+i I Zk 

\k=i+l 



fc=s+l 



- V- 

D(i,t,q i+1 , t ) 

A 3 (s,t;t)A(s,t;s+l) 



D(s,t,q s+ht ) 

split(l+, . . . , S q, (s + I)", . . . t+ (t + 1)+ . . . , n + - P~) 

s— 1 n 



(6.35) 



s/z^zn nr=i 1 (*» * + x ) 



X 



(t+ H)(ss + l) 3 f ^ _ 



i=0 j=t+l v ' J ' ^ +i >J 



\k=i+l 



- E f E - 



j=t+i 



\k=s+l 



A( g ,t;t)A 3 ( S; t; S + l) 3 



(6.36) 



7. Selected results for triple collinear limits 

To illustrate our general results for multi-collinear limits, in this section we list some 
of the triple-collinear splitting functions. The AM = splitting amplitudes are 
obtained directly from MHV amplitudes and we do not list them here. Explicit results 
are given in Section 6. For the AM = 1 (and therefore AP = 1) amplitudes, there 
are two types of splitting function corresponding to Split + (mi) and Split_(mi, to 2 ). 
In the specific case of three collinear particles, these are related by parity thereby 
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reducing the number of independent amplitudes to at most two for each splitting. 
Here, we list only the most compact form of the amplitudes. 



7.1 qgg -> q 



There are only two independent AM = 1 splitting amplitudes, which can be obtained 
by setting m = n = 3 in Eqs. (6.13) and (6.14). Explicitly we find, 



split(l+2+,3-^P+) = 



(2 3) z 2 



3/2 



split(l-,2 + ,3--.P (? -) = 



^2,3 ((2 1) (3 1) y/zi) (Z 2 + Z 3 ) 

(13)((13) V ^T+(23) V ^) 2 
"(12) (2 3) si )3 «21) v ^+(31) v ^)' 
z x (2 3) z 2 3 / 2 

V^S2,3 ((2 1) + (3 1) y/zi) (Z 2 + Z 3 ) 

as)! 

(12) (2 3) s 1)3 «21) v ^+(31) v ^)" 



(7.1) 



(7.2) 



All others can be obtained by parity and charge conjugation. These expressions 
numerically agree with the splitting functions given in [7] 4 , 



split(1 + ,2 + ,3-^P+) 



1 



x 



s l, 2^2,3 

[1 2] ((1 3) y/IT + (2 3) y/i^qi 2] + [3 2] 

S l,3 

T 2 { Zl + z 2 ) [1 2] ((1 3) y/zi + (2 3) y/zfi 



z \Z 2 Si 2 



split(l-,2 + ,3-^P,-) 



{z 2 + z 3 ) 

1 



(7.3) 



x 



S 1,2 S 2,3 

[1 3] ((1 2) ^zT + (3 2) y^) 2 ([l 3] + [2 3] 

Sl,3 

y/to( Zl + z 2 ) [1 3] ((1 2) yfz\ + (3 2) ^) 



+^^ + v^(2 3)[13] 



U2 + z 3 ) 



(7.4) 



We see that the two sets of results have the same types of singularity structure as 
z 3 — > and Zi — > 1 corresponding to the soft and double soft gluon limits. 



4 Note that there is a small typographical error in Eq. (5.56) of Rcf. [7]. s 2 z should be replaced 
by s 12 in the last term of the equation for split I"* 999 (k^ , , ) 
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7.2 qqg -> 5 

There are again only two independent AM = 1 amplitudes, both of which can be 
obtained from Eq. (6.16) by setting s = 2 and A = ±|. We find, 



split(l+,2+3-^P + ) = 



((1 2) ^il + (3 2) y/zi) ((1 3) VST + (2 3) 

(7.5) 



(12) (2 3>si )3 ((21) v ^+(31) v ^) 
z 2 (2 3) 



split(l+2- 3+-P+) = 



z7s 2 , 3 ((2 1) ^ + (3 1) ^) (23 + 2 3 ) ' 
((12) V ^T+(32) V ^) 3 
(12) (2 3)si, 3 ((21> v ^+(31) v ^) 
* 3 (2 3) 

57*2,3 ((2 1) Jzi + (3 1) (z 2 + z 3 ) ' 



(7.6) 



All others can be obtained via parity and charge conjugation. Eq. (7.5) numerically 
agrees with the analogous expression given in Ref . [7] . We were not able to find agree- 
ment between Eq. (7.6) and the expression for sp\it 9 S* 9qq (kf , k^, k£) in Eq. (5.53) of 
Ref. [7]. 

7.3 qQQ — > q and qqQ — > Q 

In this special case both colour structures lead to the same splitting amplitude. There 
is only one independent AM = 1 helicity configurations. It is given in Eq. (5.7). 



8. Conclusion 

In this paper we have considered the collinear limit of multi-parton QCD ampli- 
tudes at tree level. We have used the new MHV rules for constructing colour ordered 
amplitudes from MHV vertices to generalise our previous results for gluon-only split- 
ting functions. Our main results are general formulae for timelike splitting functions 
involving up to two negative helicity partons and an arbitrary number of positive 
helicity partons. We anticipate that the expressions presented here will be useful in 
developing higher order perturbative predictions for observable quantities, such as 
jet cross sections at the LHC or in examining the high energy limit of QCD. 

A key point of our approach is that in the collinear limit only a subset of MHV 
rules diagrams contribute - those where every propagator invariant s goes on-shell 
in the multi-collinear limit. We observe that the splitting functions have a simple 
structure, and can be written as sums over the corresponding poles in s multiplied 
by a coefficient that is either entirely composed of holomorphic spinor products (ij) 
or entirely composed of anti- holomorphic spinor products [ij]. This implies that the 
coefficients are supported on a single degree-one curve in (anti)-twistor space. 
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